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Abstract. We investigate various properties of the sublevel set {x : 
ff( x ) < 1} an d the integration of h on this sublevel set when g and h 
are positively homogeneous functions. For instance, the latter integral 
reduces to integrating hexp(-g) on the whole space R n (a non Gaussian 
integral) and when g is a polynomial, then the volume of the sublevel 
set is a convex function of the coefficients of g. In fact, whenever h is 
nonnegative, the functional J 0(g(x))/i(x)dx is a convex function of g 
for a large class of functions <f> : R+ —¥ R. We also provide a numer- 
ical approximation scheme to compute the volume or integrate h (or, 
equivalently to approximate the associated non Gaussian integral). We 
also show that finding the sublevel set {x : g{x) < 1} of minimum vol- 
ume that contains some given subset K is a (hard) convex optimization 
problem for which we also propose two convergent numerical schemes. 
Finally, we provide a Gaussian-like property of non Gaussian integrals 
for homogeneous polynomials that are sums of squares and critical points 
of a specific function. 



1. Introduction 

Positively homogeneous functions (PHF) form a wide class of functions 
that one encounters in many applications. As a consequence of homogeneity, 
they enjoy very particular properties, and among them the celebrated and 
very useful Euler's identity which allows to deduce additional properties of 
PHFs in various contexts. One goal of this paper is to show that another 
not well-known property of PHFs yields surprising and unexpected results, 
some of them already known in particular cases. Namely, we are concerned 
with PHFs, their sublevel sets and in particular, the integral 

(1.1) y^I g>h (y) := / /i(x)dx, 

•^{x:9(x)<2/} 

as a function I g ^ : M + — > R when g, h are PHFs. With y fixed, we are 
also interested in Ig t h(y) now as a function of g, especially when g is a 
nonnegative homogeneous polynomial. Nonnegative homogeneous polyno- 
mials are particularly interesting as they can be used to approximate norms; 
see e.g. Barvinok [U Theorem 3.4]. Interestingly, the integral (jl.ip is 
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related in a simple and remarkable manner to the non-Gaussian integral 
Jjg n hexp(—g)dx and therefore, any information on either one can help un- 
derstanding and evaluating the other. 

Functional integrals appear frequently in quantum Physics. In the impor- 
tant particular case when h = 1, Morosov et Shakirov pT] proved that for 
all forms g of degree k, 

(1.2) J 9il (l) = cte(fc) • / exp(- 5 )dx, 

where the constant depends only on k. In fact, a formula of exactly the same 
flavor was already known for convex sets, and was the initial motivation and 
starting point of this paper. Namely, if C C W 1 is convex, its support 
function o~c is a PHF of degree 1, and the polar C° C M. n of C is the convex 
set {x : ctc( x ) < 1}- Then 

vol(C°) = —, [ exp(-<r c (x)) dx, VC. 

Ill jRn 

(See e.g. Hiriart-Urruty [3j Exercise 65, p. 243-244] and Barvinok [4, Prob- 
lem 4, p. 207].) In fact, as we will see, the intriguing result (ll.2p is a 
particular case of a more general and striking result about PHFs. 

The connection between I g \ and J exp(-g) is a consequence of what the 
authors in |17| call action-independence of integral discriminants. That is, 
for a very general class of univariate functions (j), and for every form g of 
degree d, 

(1.3) J 0(<?(x))dx = C(cf>,d)-e(g), 

for some function where the constant C(4>,d) depends only on eft and d; 
see [T71 (73)]. And so for instance, with the choices t h-> 0i(i) = l[o,i] (*) an d 
02 (i) = ex p( — t), an d with M n as integration contour, 

f C(^,d) f 

J{x: 9 (x)<l} ^(,<P2,«j Jr™ 

whenever any of the above integrals finite. In Morosov and Shakirov [T7] the 
motivation of the authors was not to connect l g \ with J exp(-g) but rather 
to study a specific class of non-Gaussian integrals, the so-called integral 
discriminants. In fact, the theory of integral discriminants which lies in the 
field of non linear algebra, extends to integrals on more general contours of 
integration, and Ward identities permits to study properties of such integrals 
that are invariant under a change of contour of integration. Hence the goal 
in |17| is to provide exact formulas for integral discriminants in terms of 
algebraic invariants of the form g involved in the exponent; several highly 
non trivial examples are provided in [181 [22] and some of their results are 
further analyzed in Fuji [7] and Stoyanovsky [23] . 
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Contribution. In the present paper we do not study non Gaussian inte- 
grals from the point of view of algebraic invariants as in e.g. [17} [I"8l I22j . 
We are interested in properties of the integral (11.11) (and variants as well) 
and its associated non-Gaussian integral J" Rn hexp(— g)dx for PHFs that are 
not necessarily polynomials. Among other things, we are also interested in 
methods for their evaluation (or approximation). More precisely: 

(a) We first extend the action-independence property fll .3f) to PHFs of 
degree ^ d S R. Namely, given a measurable mapping <f> : R + — > R, and 
g, h PHFs of respective degree 0/d,pGl and such that J \ h\ exp(—g)dx 
is finite, 

(1.4) / ^(x))/t(x)dx = C(<f>,d,p)- f hexp(-g)dx, 

where the constant C(4>,d,p) depends only on (p,d,p. In particular, if the 
sublevel set {x : g(x) < 1} is bounded, then for every nonnegative y 6 R, 

r y{n+p)/d r 

(1.5) / hdx = — ■ — — - / hexp(-g)dx, 
J{x ■ 9(x)<j/} r(l + (n + p)/d) J R n 

with r being the standard Gamma function^. And so the (Lebesgue) volume 
of the sublevel set {x : g(x) < y} is given by: 

n/d r 

(1.6) vol ({x : 5(x) < y}) = / exp(-g) dx. 

As already mentioned, when h(x) = 1 for all x, and g is a positive definite 
form, (jl.6p was already proved in Morosov et Shakirov \17\ (75) p. 39]. In 
fact, the same arguments as in |17| can be adapted to PHFs g and with 
h ^ 1. The same technique applies for the integral 



I 9l ,...,g m ,h ■= / hdx, 

J{x:g k (x)<l, k=l,...,m} 

for a finite family (gk), h of PHFs with same degree < d G 
In addition, we obtain the identity 



{x:g(x)<j/} 



exp(- 5 ) h dx f exp(-z)z( n+p )/^ 1 dz 
_ Jo 



eM~9)hdx n(n + p)/d) 

which expresses how fast /u({x : ^(x) < y}) converges to //(R n )(= J exp(-g)hdx) 
when y — > oo and [i has the non Gaussian density hexp(-g) with respect 
to the Lebesgue measure. It is the same speed with which the truncated 
(Gamma) integral Jq ex~p(—z)z( n+p ^ d ~ 1 dz converges to T((n +p)/d). 

(b) We also provide an alternative and simple proof based on 
Laplace transform (in the spirit of Lasserre and Zeron [13] to provide 



1 The Gamma function r : R ->■ 1 is denned byoh> T(a) := f 1 " 1 exp(-t)dt, for 
every a > — 1. 
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explicit expressions for certain integrals on specific domains like e.g., a sim- 
plex or an ellipsoid). This proof permits to interpret (|1.2p as duality result in 
the usual (•, +)-algebra which parallels duality in convex optimization, i.e., 
in the (max, +)-algebra. This complements the list of the many remarkable 
properties of homogeneous functions; for instance in |12j for integration, in 
[9] for optimization, and in |14j for several properties of some conjugates of 
homogeneous functions. 

(c) Convexity: (|1.5p also helps analyze properties of the mapping fh : 
C d ->R, 

9 ^ fh{g) ■= / h(x) dx, \/ g £ C d , 

J{x:g(x)<l} 

where Cd is the space of nonnegative PHFs of degree ^ d G R with 
bounded sublevel set {x : g(x) < 1}, and h is a given PHF. Whenever h is 
nonnegative, fh is strictly convex and when g is a homogeneous polynomial, 
we provide an explicit expression of its gradient and Hessian on the interior 
of its domain. This is the analogue in our context of a formula already 
provided in Mtiller et al. |19j when the integration domain is a polytope. 

This convexity property is important. For instance, given K C M n , it 
shows that the problem of computing a homogeneous polynomial g such that 
{x : g(x) < 1} contains K and has minimum volume, is a finite-dimensional 
convex optimization problem (but hard to solve even though it is convex)! 
This problem has received a lot of attention with an elegant solution in the 
case d = 2 (that is, finding the ellipsoid of minimum volume). Notice that 
here with d > 2 the problem is still convex even though the sublevel set 
{x : g(x) < 1} is not required to be convex. We present two numerical 
approximation schemes based on semidefinite programming, that provide 
sequences of upper and lower bounds that converge to minimum volume. 

(d) Polarity: When g is a convex PHF of degree d £ M. with compact 
sublevel set G = {x : g(x) < 1}, the polar G° is a certain sublevel set of the 
Legendre-Fenchel conjugate g* of g, which is itself a PHF of degree q with 
l/d+ 1/q = 1 and therefore we also obtain the volume of G° in the form of 
non Gaussian integral, the dual analogue of the one obtained for G. 

(e) Numerical approximation: Evaluating a non Gaussian integral 
J exp(-g) is a very hard problem even if in some cases exact results can be 
obtained from some algebraic invariants of g, as in e.g. |1 T |. 1 18| . [22] . However, 
it turns out that (jl.5j) is indeed very useful to approximate non Gaussian 
integrals from polynomials. Indeed, when h is a polynomial and g is a homo- 
geneous polynomial, we show how to approximate the non Gaussian integral 
J hexp(—g)dx, as closely as desired, by solving a hierarchy of so-called semi- 
definite programs. This numerical approximation scheme complements the 
exact results mentioned above. 

(f) A Gaussian-like property: Finally we also prove a (again intrigu- 
ing) Gaussian-like property of homogeneous polynomials that are sums of 
squares (SOS). Recall that every homogeneous and SOS polynomial g of 
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degree 2d can be written as v^(x) Svrf(x) for some positive semidefinite 
matrix £ (not unique in general) and where v d (x) is the vector of all mono- 
mials of degree d. If we now consider the function 

£ ' V 0(E) = (det S) fc J exp(-fcv d (x) T £v d (x)), 

where k = n/(2d£(d)) (with t{d) = ( n+ d -1 )), then its critical points are the 
positive semidefinite matrices X) such that X -1 is the matrix of moments of 
order 2d associated with the non Gaussian density exp(— fcvrf(x) Sv^(x)) 
(normalized to be a probability density) ... exactly like in the Gaussian case 
(i.e. when d = 1) where £ _1 is a covariance matrix! (In the latter case, the 
function 6 is constant, hence all points £ are critical). 



2. SUBLEVEL SETS AND NON GAUSSIAN INTEGRALS 

2.1. Notation and definitions. Let R[x] be the ring of polynomials in the 
variables x = (xi, . . . , x n ) and let Rfxjd be the vector space of polynomials 
of degree at most d (whose dimension is s(d) := (™^ d ))- For every d G N, 
let := {a G W l : \a\ (= £i<*i) = d}, and let v d (x) = (x Q ), a G W l , be 
the vector of monomials of the canonical basis (x a ) of R[x]<2. Denote by Sk 
the space of k x k real symmetric matrices with scalar product (B, C) = 
trace (BC); also, the notation B y (resp. B >- 0) stands for B is positive 
semidefinite (resp. positive definite). 
A polynomial / G M[x]^ is written 

x->/(x) = Yl f'^i 

for some vector of coefficients f = (f a ) G M s ( d ). 

A real-valued continuous function / : R n — > R is homogeneous (resp. 
positively homogeneous) of degree k (k G R) if /(Ax) = A fc /(x) for all 
A / (resp. for all A > 0) and all x G R. For instance if x i-> /(x) is 
homogeneous of degree k, then |/| is positively homogeneous of degree k 
(but not homogeneous). Let us denote by: 

• Cd the convex cone of PHFs g of degree d G R with bounded sub- 
level set G := {x : g(x) < 1}. Notice that necessarily d ^ and 
any g G Cd is nonnegative. Indeed, if d = then x G G implies 
Ax G G for every A > 0, in contradiction with boundedness of G. 
Similarly, assume that g(xo) < for some some xo (hence xo G G). 
As g(Axo) = A d <7(xo) < for every A > 0, Axo G G for every A > 0, 
again in contradiction with boundedness of G. 

• P[x]d C R[x]d, d G N, the convex cone of homogeneous polynomials 
of degree d with compact sublevel set {x : g(x) < 1}. Hence, 
P[x] d C Cd whenever d G N. 
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Remember that a convex function is proper if /(x) > — oo for all x and 
/(x) < +00 for some x. A proper convex function is closed if it is lower 
semicontinuous; see Rockafellar [21J. 

In the sequel we will need the following result which extends to PHFs of 
degree del, one already used in [17] for forms of degree d G N (and also 
used in Hiriart-Urruty [3j for the volume of the polar of a convex set). 

Lemma 2.1. Let rf,p£R with d ^ 0, and let <j) : R+ — > R+ be a measurable 
function such that 

poo 

/ t n+p - l (i){t d ) dt < 00, 

JO 

and let g, h be nonnegative PHFs of degree ^ d £ R and p£l respectively. 
The functional / </>(g(x))/i(x) dx is /miie if and only if / foexp(— #)dx < 
+00, in which case: 

POO 

d / t n+p - eft 

(2.1) / 0(g(x))fe(x)dx = ^ • / ^exp(- 5 )dx. 
Jr™ r((n + p)/d) y R n 

V V ' 

cte(n,p,d,<p) 

Proof. With z = (z\, . . . , z n —i), do the change of variable x\ = t, xi = 
tzi, . . . , x n = tz n -i) so that one may decompose L n <j>(g(x))h(x)dx into the 
sum 

f t n+p - 1 ^{t d g(l,z))h(l,z)dtdz 
+ / t n+p - 1 4>(t d g(-l,-z))h(-l,-z)dtdz, 

t n+p - x (t){t d g{l,2))dt \ h(l,z)dz 
+ f ( t n+p -^(t d g(-l, -z)) dt] h(-l, -z) dz, 

JR"- 1 \J0 J 

where the last two integrals are obtained from the sum of the previous two 
by using Tonelli's Theorem (since <f> and h are nonnegative functions) and 
whether or not the integrals are finite; see e.g. Dunford and Schwartz [6j 
Theorem 14, p. 194]. 

Next, with the change of variable u = t g(l, z) l l d and u = tg(—l, —z) l l d 

(2.2) / </>(g(x))h(x)dx = ( f u n+p - 1 ( f ) (u d )du).A(g,h), 
Jr" \Jr + / 

with 

(2-3) A(g, h)=! ( , H l; Z } )/H + f h{ -\: Z } VH ) dz. 
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In particular, the choice t \-t 4>(t) = exp(— t) yields: 

d f 

(2.4) A(g,h) = / 7iexp(-p)dx, 

and so A(g,h) is finite if and ony if J hexp(—g)dx is finite. Substituting 
A(g, h) in $T2$ with yields (j2TTj) . □ 

2.2. Integration on sublevel sets and non Gaussian integrals. As a 

consequence of Lemma 12. II we obtain the following result. 

Theorem 2.2. Let g, h be PHFs of degree ^ d £ R aniip G R respectively, 
and such that g £ C^. 27ien /or every y G [0, oo): 

(2.5) vol({x : p(x) < y}) = — . . / exp(-p) dx, 

r(l + n/d) J R n 

and 

r y( n +p)/ d r 

(2.6) / hdx = =— -, r-r— I /i exp(— a) dx, 

•/{*:*(*)<«} r(l + (7l+p)/d) An PV ^ 

whenever J Rn exp(— p)dx (or J/ X . g ( x )<n |/i|dx,) is finite. And 112. 6}) also 
holds if h is replaced with \h\, or with max[0, h], or with max[0, — h]. 

Proof. Write h = h + — h~ with h + (x) : = max[0, h(x)], h~(x) := max[0, —h(x)] 
for all x. Observe that \h\, h + and h~ are nonnegative and positively ho- 
mogeneous of degree p, and \h\ = h + + h~ , h = h + — h~ . 

With A as in ()2.2h . using Lemma 12 . 1 1 with (f)(t) = exp(— t), yields 

in / s , r (l + (n +p)/d) .. ., .. 

/j exp -p dx = v • A(g, \h\) 

n + p 

where all integrals are finite since the right-hand-side is finite. Similarly, 
using Lemma ET] with 4>{t) = l[o,i](i), yields 



f \ h \dx = -J—.A(g,\h\), 

•/{x : g(x)<l} n +P 



and therefore, 



/ \ h \ dx = r/i i I 1 i \ia\ / I' 1 ' ex P(-9)dx. 

iix : g(x)<n r(i + in + p)/d) y R n 



'{x:g(x)<l} 

Next, by homogeneity, 

f \h\dx = y( n +Py d [ \h\dx, 

•/{x: 9 (x)<s/} /{x: ff (x)<l} 

and so 

(2.7) / \h\dx = y / \h\ exp(-p)dx. 

J{x : 9 (x)<y} I\l + (n + p) / d) J R n 
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Obviously, with same arguments, (12. 7h also holds if we replace \h\ with h + 
and h~ and finiteness follows because 

< / h + dx < \h\dx, 

and similarly for h~ . Hence (|2,7p with h in lieu of \h\ (i.e. (|2.6p ) follows by 
additivity since h = h + — h~ . 

Finally, from what precedes, finiteness of Jr x . fl ( x )<n \h\dx is equivalent to 
finiteness of L n \h\ exp(— g)dx. □ 

When y = 1, g is a positive definite form with d £ N, and /i is identical 
to the constant function 1, (|2.5p is already proved in Morosov and Shakirov 
13 (75) p. 39]. 

In particular, when g is the quadratic form x i— > Q2 (x) := ^x r Qx for some 
real symmetric positive definite matrix Q, one retrieves that the volume of 
the ellipsoid £(?/) := {x : 52 (x) < y} is simply related to the determinant 
of Q by the formula 

(2.8) vol(g(y)) = * / exp(--x T Qx)dx 

r(i + n/2) y R n 2 

V v ' 

=(2 7 r)"/ 2 /v / dStQ 

So Theorem 12.21 states that the volume of the sublevel set is simply related 
to the integral of exp(— g(x)) over the entire domain R n which happens 
to be simply related to the determinant of Q when g is the quadratic form 
x r Qx. One goal of the theory of integral discriminants is precisely to express 
J exp(-g) in terms of invariants of g when g is a form. See e.g. Morosov 
and Shakirov |17^ I18j and Shakirov |22j . 

Intersection of sublevel sets. Suppose that with h being positively ho- 
mogeneous of degree p, one wishes to compute the integral J fi h(x) dx where 

fi:={x£l" : g k (x) < z k , k = l,...,m}, 

for some m PHFs g±, ... , g m of degree 7^ d € K, and some (strictly) positive 
vector z £ M m . Equivalently, 

Q:={x£ R n : g k (x; z) < 1, k = 1, . . . , m}, 

for the functions x 1— > gk{x; z) := gu{Zu l ^ d x), k = 1, . . . , m, which are also 
PHFs of degree d G IR. Hence with no loss of generality, one may restrict to 
sets of the form 

(2.9) Q(y) := {x : g k {x) < y, k = 1, . . . , m} 

for some positive scalar y S K, and PHFs gi, . . . , g m of same degree dsK. 
Notice that x h-> ^(x) := max bi( x )) • • • ; 9m (x)] is a PHF of degree d. 
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Corollary 2.3. Let h be a PHF of degree let g%, . . . , g m be PHFs of 

degree ^ d S R, and assume that the set {x : <?fc(x) < 1, fc = 1, . . . , m} is 
bounded and for every y > 0, Ze£ f2(y) fre as in $2.9\) . Then: 

(2.10) vol(fifo)) = T( f' d M , I exp(-V)dx, 

r(i + n/d) 7 K n 

and 

/• y( n +P)/ d r 

(2.11) / /i(x)dx = / ft exp(-V>) dx, 

yno,) r(i + (n + p)/d) y R n 

whenever J" Kn exp(— -0) dx. 

Proof. Notice that x 1— >■ V( x ) := max [<?i( x )) • • • > 9m ( x )] i s a l so a PHF of 
degree / d G 1, and O(y) = {x : V( x ) — ?/}• And so if is bounded 
then V G C d . Hence (|23IJ) - (|ZI0] 1 is just ([23 ]> -(|23 ]) with h and V in lieu of 
h and 5. □ 

2.3. An alternative proof with a duality interpretation. Next, we 
present an alternative proof of Theorem 12.21 that uses Laplace transform 
techniques and provides an interpretation of the result in an appropriate 
duality framework. 

Suppose that g,h G C^. Since g is nonnegative, the function I g h vanishes 
on (—oo,0]. Its Laplace transform Cj g h : C — > R is the function 



roc 

A ^ C I 9 ,h( X ) : = I eX P 
JO 



(-><y)ig,h(y) dy, 



and observe that 



£ighW = / ex P(-^) / hdic) dy 

JO \J{x:g(x.)<y} J 

= / ft(x) / exp(— Xy)dy ] dx [by Fubini's Theorem] 
7m« \d ff ( x ) / 



1 

A 

X-P/d 



A 



ft(x) exp(— A<?(x)) dx 

h(X 1 ^ d x) exp(-5(A 1 / d x)) dx [by homogeneity] 
l_ w ^/ l ( z ) exp (- 5 (z))d Z [b y A-^ x ^z] 



A!+( 

ft(z) exp(— g(z)) dz 



r(l + (n + p)/d) 
And so, by uniqueness of the Laplace transform, 

(n+p)/d 



y(n+p)/a r 
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which is the desired result. So the above expression of I g< h(y) is obtained 
by "inverting" the Laplace transform Ci g h at the point y, which in fact, as 
we next see, is solving a "dual" problem. 

For analogy purposes, consider, the optimization problem 

Pg,h{y) = sup{/i(x) : y(x) < y}, 

x 

where y is fixed, h and —g are concave and h is nonnegative. Equivalently, 
Pg,h{y) = exp(0 fli 7j(?/)) where y h-> 6 9t h(y) is the optimal value function of the 
optimization problem 

: 9g,h(v) = sup{ln/i(x) : y(x) < y}. 

X 

Associated with P y is the cfata/ problem P* : mf\{Q(X) : A > 0}, where 
G : R + ->• R is the function A ^ G(A) := sup x {ln/i(x) + A(y - y(x))}. 
Observe that: 

A (->■ G(A) = Ay + sup{ln/i(x) - Ay(x)} 

X 

= Ay + sup { ln(A- p / fl! /i(z)) - y(z)} [via z = A 1 / d x] 

z 

= Ay H — — In A + sup{ln h(z) — y(z)}. 

Ct z 



And so the dual problem P* reads 

4- 

a>o d 



P* : 7 = sup{ln/i(z) - y(z)} + inf {Ay - ^ In A} 



= lny p / d + In ^sup{/i(z) exp(-y(z))}^ + |(1 - In 



In particular, if /i is log-concave and g is convex, by a standard argument of 
convex optimization, 7 = 9g t h(y) (= ln/o fli ^(y)). And therefore, 

p,, fc (y) = exp(VG/)) = ^^f^T sup{/ l (x)exp(-y(x))} 

to compare with 

I,M=V^« WT ^ Tm I Mx)exp(- 9 (x) )(J x. 

Alternatively, the Legendre-Fenchel transform of g ,h(y) (fo r the concave 
version) is the function 

0g,hW = mf{Ay-V(y)} 

= inf{Ay - sup{ln h(x) : y(x) < y} } 

= inf {Ay(x) — ln/i(x)} = — In I sup{/i(x) exp(— Ay(x))} 
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and so when h is log-concave and g is convex, 6g t h(y) = (d* /l )*(y), so that 
g , h (y) = m{{\y-e: h (X)} 



inf {Ay + In ( sup{h(x) exp(— Ag(x))} ) } 



= ln^/ d + lnr S up{^( Z )exp(- 5 ( Z ))}") + |(l-ln|). 
In summary, to the Laplace transform step 
(2-12) C Igh (X) = [ exp(-\y)I g;h (y)dy 

in the usual (•, +)-algebra, corresponds to the Legendre Fenchel transform 

°* g ,hW = inf {-% - OgMy)} = - su p{- A y + QgAv)} 

y y 

= -lnsup{exp(-Ay)p 9i/l (y)}, 

y 

where to emphasize the analogy, the latter term can be written 

lnf" [ n exp(-><y) Pg,h(y) 



i.e., the "sup" operator is integration "J" in the (max, +)-algebra. And 
with this convention 



exp(6»* h (A)) ~ "J n exp(-Xy) Pg , h (y). 
(Compare with (|2.12p .) Similarly, the Laplace inverse transform step 

(2-13) Ig,h(y)= exp(Xy)£ Igh (X)dX 

J ui—ioo 

(where ui G C is on the right of all singularities of Ci h ) and which yields 

y(n+p)/d r 

r(l + (n + p)/d) J Rn 

is solving the "dual" and corresponds to the Legendre-Fenchel transform 
(which is involutive) applied to 8* h 

exp(V(2/)) = exp(inf{Ay - 0* h (A)}) ~ -" / '" exp(-Ay) exp(^ ft (A)) dX 
(compare with (|2.13p ). and which yields 

eMOgM) = y P/d ^07 "/" Mx)exp(- 5 (x)). 

A rigorous analysis of the links between integration and linear optimization 
on a polytope has been already investigated in [ID] . 
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2.4. Approximating non Gaussian integrals. As we have already men- 
tioned, in some cases the non Gaussian integral can be computed explicitly 
in terms of some on algebraic invariants of g; see e.g. Morosov and Shakirov 
[17j and Shakirov [22] . But so far there is no general formula and therefore 
an alternative is to seak for a numerical scheme for its evaluation, or at least, 
its approximation. 

For this purpose, we next show that Theorem 12.21 is helpful as it provides 
a means to compute any moment of the measure dfx = exp(—g)dx on W n 
by computing the same moment but now of the Lebesgue measure on the 
sublevel set {x : g(pc) < 1}. Indeed, for every a E N n , letting x i— > h(x) := 
x a in Theorem 12.2] yields 

(2.14) f x a exp(-g(x))dx = r(l + (n+|a|)/d) / x a dx, 

JB" J{x:g(x)<l} 

where \a\ = ^i a «- Have we made any progress with this equivalence? 

The answer is yes. If g is a (non necessarily homogeneous) polynomial, 
it turns out that every moment of the Lebesgue measure on the sublevel 
set {x : g(x) < 1} can be approximated as closely as desired by solving a 
hierarchy of semidefinite programs^ as described in Henrion et al. [8|. In 
fact, for every a E N™ fixed, the moment 

z a := / x a dx 

•/{x: 9 (x)<l} 

can be approximated to arbitrary precision e > fixed in advance, by solving 
two sequences of semidefinite programs, one which provides a monotone non 
decreasing sequence of upper bounds Uk, k E N, while the other provides a 
monotone non increasing sequence of lower bounds £k, k E N. The procedure 
stops whenever ut — Ik < e, in which case one may set 

z a w z a := (u k +4)/2- 

This requires to solve two sequences of semidefinite programs for each a E 
N n . In fact, if one is ready to relax the monotonicity property of the upper 
and lower bounds {uk,lj~}, it is enough to solve a single sequence of semi- 
definite programs, e.g., the one defined to approximate the mass Zq. Then 
if d E N and e > are fixed, and k is large enough, not only \uk — zq\ < e 
but also from the solution of the semidefinite program at step k one obtains 
scalars z a such that \z a — z a \ < e, for all a E N n such that \a\ < d. How- 
ever, in contrast to the case of upper and lower bounds, there is no simple 
stopping criterion to guarantee the e-approximation. For more details, the 
interested reader is referred to Henrion et al. [8J. And therefore, once the 



A semidefinite program is a finite-dimensional convex conic optimization problem, 
that up to arbitrary (fixed) precision, can be solved efficiently, i.e., in time polynomial in 
the input size of the problem. For more details the interested reader is referred to e.g. 
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have been computed, using Theorem 12.21 we obtain: 



/ x a exp(-o(x))dx - z a T(l + (n+\a\)/d) 



< e, Va G N% 



which provides an e-approximation guarantee for the non Gaussian integral 
L„ exp(— g)dx (and more generally for the integral J Rn hexp(—g)dx when- 
ever h is any polynomial). 

2.5. Sensitivity analysis and convexity. Recall that when d G N, P[x]d C 
C d is the convex cone of nonnegative and homogeneous polynomials of de- 
gree d, with compact sublevel set {x : g(x) < 1}. Formula (|2.6p of Theorem 
12.21 allows us to provide insights into the function / : C d — > R, defined by: 

(2.15) 9^fh(9)-=f /i(x)dx, g£C d , 

•/{x:g(x)<l} 

where /i is a PHF. In particular when one wishes to see how fh changes 
when some coefficient of g G P[x]d varies. Notice that the restriction of fh 
to P[x]d may be seen as a function fh : R^) — > R of the coefficient vector 
of the polynomial g G P[x]<2, where £(d) = ( n+ d ~ 1 )- 

Before proceeding further we need the following result. Recall that the 
support supp /i of a Borel measure on R n is the smallest closed set A such 
that /i(R n \ A) = 0. Let C° := {g G C d ■ g is continuous} C C d . 

Lemma 2.4. Let /i be a non trivial a-finite Borel measure on R n and let 
G^ C R[x]d be the convex cone of polynomials g of degree at most d such 
that f exp(— g)d/i < oo. Then: 

(a) With d G R, the function f : C d — > R, with g h-» /(g) := j exp(—g)dfj,, 
is convex (and strictly convex on C° z/supp/i = W 1 ). 

(b) If /i(0) > /or some open set O C R n and if d G N, £/ien / is strictly 
convex and twice differentiable on int(Q^), with: 

(2.16) = y x a exp(- 5 )d/i, Va, |a| < d 

(2.17) = / ^ exp(- ff ) d//, Va,/3, |a|, |/3| < d. 
og a ogp J 

Proof, (a) Observe that / is nonnegative. Let a G [0, 1] and let g, q G C^. 
To prove /(ag + (1 — a)q) < af(g) + (1 — a)f(q), we only need consider the 
case where f(g), f(q) < +oo, for which we have 

f(ag + (1 - a)q) = j exp(-ag - (1 - a)q) dfi. 

By convexity of u \-t exp(— u), 

f(ag + (l-a)q) < j [a exp(-g) + (1 - a) exp(-g) ] dfj, 
= a f(g) + (l-a)f(q), 
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and so / is convex. Now, in view of the strict convexity of u t— >■ exp(— it), 
equality may occur only if g(x) = g(x), ^-almost everywhere. If g, q G C^, 
the set A := {x : g(x) — q(x) = 0} is closed and so if A / M n then 
/i(A) < ^(M n ) because supp/i = W 1 . Therefore, equality occurs only if 
g = q so that / is strictly convex on C%. 

(b) Next, if d G N and g G int(0u), write g in the canonical basis as 
f( x ) = E Q 5« xa - F ° r every a G NJ}, let e Q = (e a (/3)) G be such that 

= 5/3= Q (with 5 being the Kronecker symbol). Then for every t > 0, 



Notice that for every x, by convexity of the function t \— > exp(— ix a ), 



because for every x, the function t \— > ip(t,x) is nondecreasing; see e.g. 
Rockafellar |21|, Theorem 23.1]. Hence, the one-sided directional derivative 
f'(g;e a ) in the direction e a satisfies 



where the third equality follows from the Extended Monotone Convergence 
Theorem [21 1.6.7]. Indeed for all t < t$ with to sufficiently small, the 
function if>(t, •) is bounded above by tp(to, •) and J exp(—g)ip(to,x)dfi < oo. 
Similarly, for every t > 




ip(t,x) = inf^(t,x) = exp(— tx Q )j 4=0 



= — x 



.a 





and by convexity of the function 1 *— > exp(tx a ) 




= x' 



a 



Therefore, with exactly same arguments as before, 



\ ,. f{g -te a ) - f(g) 
f [g; -e a ) = hm 




f'(g;e a ) 



and so 
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for every a with \a\ < d, which yields (|2.16p . Similar arguments can used 
for the Hessian V 2 f(g) which yields (|2.17p . 

So the Hessian V 2 f{g) is the matrix E S^) whose rows and columns 
are indexed in the set := {a £ N n : \a\ = d} and with entries 

M d (a, P) = J x a+ ^ exp(-g) dp , a, p € r d , 

dv 

i.e., Mrf is the matrix of 2d-moments of the finite Borel measure v. Let 
h € B/W be the coefficient vector of a non trivial and arbitrary homogeneous 
polynomial h £ M[x] d . Then 

<h,M d h) (= f h{x) 2 dv{x) J > 

because n(0) > (hence v{0) > 0) on some open set Ocl" Therefore 
V 2 /(g) >- 0, which in turn implies that / is strictly convex on int(O^). □ 

Corollary 2.5. Let h be a PHF of degree p G R and wit/i / d £ R, 
consider the function fh : — )■ R defined by: 

(2.18) 5 ^ fh{g) ■= ! h(x)dx, Vg£C d . 

J{x:g(x)<l} 

T/ie function fh is a PHF of degree —(n+p)/d and convex whenever h is 
nonnegative (and strictly convex if h > onl"\{0} / ). In addition, ifd£~N 
and h is continuous with J \ h\ exp(—g)dx < oo, then: 

(a) fh is twice differentiable on int(P[x] d ) ; and for every a, ft £ N^' 

(2 - 19) ^ = ra + ^rt/eo l^^^-^^ 

(2.20) = -J; (2+ / n + ^ / x^(x)^x. 

V ' r(l + (n+ P )/d) y { x: g( x)<l} ^ ^ 

(2.21) = 1 / x-+^(x)exp(- 5 (x))dx. 
" ' ■ r(l + (n + p)/d) 7 M n 



(7>j is non trivial and nonnegative, then fh is strictly convex on int(P[x] d ) ; 
and its Hessian V 2 // l (g r ) is the matrix of 2d-moments of the measure 

u(B) = ———^— — — - / hexp(-g)dx, BeB. 
T(l + (n + p)/d) J B 

Proof. With A > and g £ C d , 



fhi^g) = / hdx= / hdx., 

i{x:Ac/(x)<l} J{x:g{X 1 / d x)<l} 

and so, doing the change of variable z = A 1 / o! x, one obtains fhi^d) 
\~ {n+p)/d f h (g), i.e., / h is a PHF of degree -{n + p)/d. 
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Next, first consider the case where h is nonnegative, and let \x be the 

cr-finite measure defined by /x(-B) := / h(x)dx for every Borel set B of R™. 

J B 

By Theorem 12.21 whenever g £ C d and fh(g) is finite, 

U9) = T { l + (l + P )/d)L^- 9) ^ 

Hence by Lemma [2.4l fh is convex (and strictly convex if h > on M n \{0}). 

In addition, if d G N and h is continuous, //1(g) < oo if g G int(P[x]d) (and 
so (7 G int(G^)). Moreover, /i being non trivial, nonnegative and continuous, 
h > on some open set O and so /J,(0) > 0. Therefore, by Lemma l2~iT b). 
fh is twice differentiate and strictly convex on int(P[x]^) and (|2.16|) - (|2.17p 
yield <^W\ and (I2T2T]) . while (f^UD follows from ([ZTgD and Theorem El 
That V 2 fh(g) is the matrix of 2a!-moments of the measure a!^ = hexp(—g)dx 
follows from (12.21 H . This proves (b). 

To prove (a) when h is not nonnegative, write h = h + — h~ with h + : = 
max[0, h] and h~ := max[0, — h]. Both h + and h~ are continuous PHFs of 
degree p, and nonnegative. Moreover, fh = f h + — f h - and so applying (b) 
to fh+ and fh-, yields (a) by additivity. 

□ 

Remark 2.6. (a) Notice that proving convexity of fh directly from its defi- 
nition 12.18]) is not obvious at all whereas it becomes much easier when using 
Theorem \2.SX 

(b ) In Lemma \2.4\ differentiability of f on the convex cone should be 
now in the sense of Gateaux- differentiablity, not explored here. 

We end up with the following relatively surprising results which even 
though are again particular cases of Lemma 12.11 deserve special mention. 

Lemma 2.7. Let y > be fixed, let h be a PHF of degree p G R and let 
£, tjj : R + — 7- R be measurable functions such that 

f t n+p -^(t d )dt < +oo. 

Let fh : Cd — > R, / d G R, be the function: 

9^fh{g)-= / C(^(x)) h(x) dx, g£C d . 

J{x:V( 9 (x))<|/} 

Then whenever f \h\ exp(—g)dx < +oo, fh(g) is finite, and 

d ! t n+ P-^(t d )dt 
(2-22) f h (g) = W)^> f heM _ g)dx , 

V v ' 
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where the constant depends only on ^,ij),d,p,y and neither on g nor h. 
Therefore, fh is convex whenever h is nonnegative (and strictly convex on 
C® whenever h > on W 1 \ {0} ). In particular, 

(2.23) / gh exp(-g) dx = — - — / h exp(-g) dx 

(2.24) / ghdx = — — \T3\ [ hexp(-g)dx 
Jfx:q( x )<i} T((n + p)/d) J Rn 



'{x:g(x)<l} 

/ exp(- 



y 

exp(-g)dx I exp(-z)z n/d ~ l dz 



T{n/d) 



\ exp(-g)dx 

/ exp(-g)dx 

(2.26) / exp(g)dx = ^— / exp(z)^ 1 dz. 

•/{x:g(x)<y} L[n/d) J 

Proof. We first assume that h is nonnegative and f hexp(-g) < +oo, so 
that <^T2\i follows from Lemma O with t H- 4>{t) := ^(t)I^ y ](tp(t)), which 
yields fh(g) = cte(y,p,d,(,,Tp) ■ A(g,h), with A(-,-) as in (|2.3p and 

cte(y jP ,d,^) = [ t n+ v~ l t(t d )dt, 

and the result follows by recalling that with <p{t) = exp(— t) one had already 
obtained in ([27 



A (9,h) = — — ■ — -— / h exp(-g)dx(< +oo). 
T((n + p)/d) 7 R „ 

When h is not nonnegative, writing h = h + — h~ where both h + and h~ are 
also PHFs of degree p, the result follows by additivity since J \ h\ exp(—g)dx < 
+oo only if both J h + exp(— g)dx and f h~ exp(—g)dx are finite. 

Finally, (I2.23p - (I2.26I) are special cases of (I2.22|) with respective choices 
t t/j(t) := 0,£(t) := exp(-t), then t \-> ip{t) := t,£(t) = t and finally, 
1 1->- rp(t) = t,£(t) := exp(-t) and t \-t £(t) := exp(t). 

At last, when h is nonnegative, convexity and strict convexity follow from 
Lemma 12.4( a) with d/i = hdx (and so supp \x = W 1 if h > 0) . □ 

So Lemma [2 .71 shows that / is convex provided that h is nonnegative and 
no matter how the functions £ and ip behave! 

Next, if n is the non Gaussian measure d[i = exp(—g)dx on R n , then 
(|2.25|) shows how fast /i({x : g(x) < y}) converges to the non Gaussian inte- 
gral f Rn exp(—g)dx as y — > oo. It converges as fast as the one-dimensional 
integral / y t n l d ~ l exp(-t)dt converges to the Gamma function T(n/d). 
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2.6. Polarity. We here investigate the polar G° set of the sublevel set G := 
{x : g(x) < 1} assumed to be compact and when g is a proper closed convex 
PHF. In this case G is a convex body, and in fact g is a gauge. 

The polar C° of a set C C R n is the convex set defined by 

C° = {x G R n : ctc(x) < 1} with <r c (x) := sup{(x,y) : y G C}. 

y 

and (C°)° is the smallest convex balanced set that contains C. 

Recall that the Legendre-Fenchel conjugate f* : R n ->RU {+oo, — oo} of 
/ : R n -> R is defined by 

/*(u) := sup{(u,x)-/(x)} uel" 

x 

The conjugate g* of a PHF 5 of degree d G R is tself a PHF of degree q with 
3 + - = 1 (where d does not need to be positive); see e.g. Lasserre [T5] . 

Proposition 2.8. Zei 5 be closed proper convex PHF of degree 1 < d G R, 
and Zei G := {x : g(x) < l/d}. Then with g + ~ = 1, 

(2.27) G° = {x G R n : #*(x) < 

where g* is the Legendre-Fenchel conjugate of g. In other words, G° is a 
sublevel set the PHF g* of degree q. Moreover, if G is bounded then 

(2.28) vol(G°) = 1 / exp(- 5 *)dx. 

q n /iT(l + n/q) J R n 

Proof (123?) is from Rockafellar [21, Corollary 15.3.2] and (I23H]) follows 
from Theorem 12.21 applied to the PHF g* and with y = 1/q. □ 

Example 1. Let x 1— > g(x) := |x| 3 if x > and +00 otherwise. Hence, 
d = 3, q = 3/2, G = [-1,1], a G (x) = \x\, and g*(x) = ^\x\ 3/2 . One 
retrieves that G° = G = [-1, 1] = {x : g*(x) < (d - l)/d q }. 

Example 2. Let x h-> g(x) := \x\ so that G = [—1,1]. As g*(x) = if 
x G [—1,1] and +00 otherwise (a PHF of degree 0) one may check that 
indeed G° = [—1,1] = G and holds although d = 1 and g is not 

strictly convex. 

Example 3. Let x \-t g(x) := xf + x\ and G = {x : x\ + x\ < 1/4}. 
Then g*(x) = 3(x\ /3 + x^ 3 )/^ 3 (a PHF of degree A/3) and G° = {x : 
^ /3 + 4 /3 <l/4V3}. 

2.7. A variational property of homogeneous polynomials. We end up 

this section with an intriguing variational property of homogeneous polyno- 
mials that are sums of squares. 

Let Vrf(x) be the vector of all monomials (x a ) of degree d and let g G 
R[x]2rf be homogeneous and a sum of squares, that is, 

(2.29) g(x) = --U d (x) T £v d (x), x£i", 
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for some real symmetric 1(d) x 1(d) matrix £ which is positive definite 
(denoted £ y 0). If d = 1 it is well-known that 

(2vr)"/ 2 

exp( - 9)dx = 7^' 

and 

/ v d (x)vj(x)exp(- 5 )cix = V r-— V-\ 
jR n vdet 2j 

that is, £ _1 is the covariance matrix associated with the Gaussian proba- 
bility density (2vr)- n / 2 v / detE exp(-g). 

When d > 1 the non Gaussian integral can be still expressed as a (possibly 
complicated) combination of several algebraic invariants of g, but in general 
not in terms of the single algebraic invariant det £. 

It turns out that det SI and f exp(— |v,i(x) T £vd(x))<ix are still related 
in the following Gaussian-like manner. Let <S^^ C S^) be the convex cone 
of £(d) x 1(d) positive definite matrices, and let 6^ : Sj/t — > R be defined by: 



(2.30) d (£) := (det£) fc f exp(-A;v d (x) T £ v d (x)) dx, 

JR n 

where k = n/(2d£(d)) and let 

/ v d (x)v d (x) T exp(-/cv d (x) T £ v d (x)) dx 



M d (E) := 



/ exp(-fcv d (x) T £ v rf (x))cix 

JR™ 



be the matrix of moments of order d, associated with the non Gaussian 
probability measure 

/ exp ( - kvd (x) T £ v rf (x) ) cix 

(2.31) n{B) := If , BeB. 

/ exp(-A;v d (x) T £v (J (x)) dx 

Observe that 8^ is nonnegative and positively homogeneous of degree 0; 
therefore 9d is constant in any fixed direction £. In particular, if d = 1 
then k = 1/2, fi is a Gaussian probability measure, Mi(£) is the associated 
covariance matrix £ _1 and #d(£) is constant. In fact, 

Lemma 2.9. Lei d 6e the function defined in $2.30\) . Then (M d (£), £) = 
1(d) for all £ in the domain of Qd an d V#^(£) = if 



(2.32) M d (£) = £ 



-i 
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Proof. Let g(x) = fcv"d(x)Svd(x) so that 



(Md(S),S) / exp(- 5 )dx = 




v<i(x)v d (x) T exp(-#) dx, 
/cv d (x) T £v d (x) exp(-g) dx 



5 exp(-s) dx 



exp(- 5 )dx [bv ([253]) ] 



which yields the desired result (M(j(S), I]) = 1(d). Next, write the gradient 
V9(E) in the form + A 2 with 



= -A;M d (£)0 d (E). 
This yields A 1 +A 2 = fc^(5])(5]" 1 -M d (S)) and so V0 d (£) = Oif M d (£) = 



Lemma 12.91 states that for all critical points XI, or equivalently for all 
critical SOS homogeneous polynomials g of the function 9^ (assuming that 
at least one such critical point exists), their associated non Gaussian measure 
djj, = exp(— g)djj, (rescaled to a probability measure) has the Gaussian-like 
property that is the "d-covariance" matrix of 

3. SUBLEVEL SET OF MINIMUM VOLUME CONTAINING A COMPACT SET 

If K C M. n is a convex body, computing the ellipsoid of minimum volume 
that contains K is a classical problem which has an optimal solution called 
the Lowner- John ellipsoid; see e.g. Barvinok [H p. 209]. In this section we 
consider the following generalization: 

P: Find a homogeneous polynomial g of degree 2d such that its sublevel 
set G := {x : p(x) < 1} contains K and has minimum volume among all 
such sublevel sets with this inclusion property. 




where E A is the adjugate of E (see e.g. [H p. 411]), and 




S" 1 . 



□ 
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With K C R n , let C2d(K) C K[x]2d be the convex cone of polynomials of 
degree at most 2d that are nonnegative on K. Recall that P[x]2d C M[x]2d 
is the convex cone of homogeneous polynomials of degree 2d with compact 
sublevel set {x : g(x) < 1}. We next show that problem P is a convex 
optimization problem: 

Proposition 3.1. The minimum volume of a sublevel set {x : g{x) < 1}, 
g G P[x]2d; that contains K C K n is p/T(l +n/2d) where p is the optimal 
value of the finite- dimensional convex optimization problem: 

(3.1) V: p= inf {/ exp(- 5 ) fix : 1 - g G C 2d (K) ) . 

Proof. From Theorem! 



vol({x : g(x) < 1}) = 1 / exp(-g)dx. 

r(l + n/2d) J R n 

Moreover, the sublevel set {x : g(x) < 1} contains K if and only if 1 — g G 
C2d(K), and so p/T(l + n/2d) in (|3.ip is the minimum value of all volumes 
of sublevels sets {x : g(x) < 1}, g G P[x]2<2, that contain K. Now since 
9 ^ Jmn exp(— g)dx is strictly convex (see Corollary 12.5( b)) and C2d(K) 
is a convex cone, problem V is a finite-dimensional convex optimization 
problem. □ 

We also have the following characterization of an optimal solution of V 
when it exists. Let M(K) be the convex cone of finite Borel measures on K. 

Theorem 3.2. Let K C M. n be compact and consider the convex optimiza- 
tion problem V in ( f3.il) . 

(a) Suppose that g* G P[x]2d is an optimal solution of V . Then there 
exists p* G M(K) such that 

(3.2) f x a exp(-g*)dx = [ x a dp* , V|a| = 2d; f dp* = 0. 

In particular, p* is supported on the real variety V := {x G K : g(x) = 1} 
and in fact, p* can be substituted with another measure v* G M(K) supported 

on at most ( n+ 2d _1 ) + 1 P°i n t s ofV. 

(b) Conversely, if g* G P[x] 2 d and p* G M(K) satisfy \3.S!\) then g* is an 
optimal solution ofV. 

Proof, (a) We may and will consider g* as an element of R[x]2d with g^ = 
whenever \/3\ < 2d. As K is compact, there exists 6 G ~P[x]2d such that 
1 — 9 G intC2d(K), i.e., Slater's condition holds for V. Indeed, choose 

■= M _1 ||x|| 2d for M > sufficiently large so that 1 - 6 > on K. Hence 
with \\g\\i denoting the ^i-norm of the coefficient vector of g (in M[x]2d), there 
exists e > such that for every h G B(6, e)(:= {h G M[x] 2 d : ||0 - h\\i < e}), 

1 - h > on K. 
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Therefore, the optimal solution g* satisfies the KKT-optimality condi- 
tions, which read: 

(3.3) / x a exp(- 5 *)dx = y*, V|a| = 2d; (l-g*,y*)=0, 

for some y* = (y*), a G N£ d , an element of the dual cone C 2 d(K)* C R S ( M ) 
of C 2d (K). By Lemma O in MM 

C 2d (KT = {y G R S W : 3f i* G M(K) s.t. y a = / x a c^*, a G N^}, 

Jk 

and so (|3.2|) is just (|3.3|) restated in terms of fi*. Finally, the last statement 
follows from Anastassiou [U Theorem 2.1.1, p. 39] applied to the ( n+ 2d~ 1 ) 
equality constraints of (|3.2p . 

(b) As Slater's condition holds for V, the KKT-optimality conditions (I3.2p 
are sufficient to ensure that g* is an optimal solution on V . □ 

Theorem 13.21 states that when ()3.2[) holds, there is an optimal solution g* 
to V such that g(x) = 1 on (at most) ( n+ 2d~ l ) + 1 points of K, the analogue 
for d > 1 of the well-known property of the Lowner-John ellipsoid in the 
case d = 1. 



Even though being convex and finite-dimensional, V is by no means easy 
to solve because there is no simple and computationally tractable way of 
describing the convex cone C2d(K). However, there are cases where one may 
provide a sequence of inner or outer approximations that both converge to 
C2d(K). One such case is when one knows all moments of a finite Borel 
measure whose support is exactly K, and another case is when K = {x : 
9fc(x) > 0, k = l,...,m} for some polynomials (5%) C R[x], i.e., K is a 
compact basic semi- algebraic set. 

3.1. Lower bounds via inner approximations. Suppose that one knows 
all moments z = (z a ), a G N n , of a finite Borel measure \i on K, i.e., 

z a = I x a d/i(x), Va G N n , 

whose support is exactly K. For every k G N and p G M[x], let Mfc(p, z) 
be the localizing matrix with respect to the polynomial p and the moment 
sequence z, that is, Mfc(p, z) is the s(k) x s(fc) real symmetric matrix with 
rows and columns indexed in the canonical basis (x a ), a G N^, of M[x]fc, 
and with entries 

M k (p, Z )(a,p) = [ p(x)x Q +^(x), Va,/3GN£ 

JK 

= ^ P7 z a+/3+j 
7 

(when p(x) = ^-7^7 x7 )- We recall the following result. 
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Lemma 3.3. ( [15} Theorem 3.2]) Let K C K n be compact and let p be a 
finite Borel measure with support K and with moments z = {z a ), ct G W 1 . 
Then p G M[x] is nonnegative on K if and only if Mfc(p, z) ^ for every 
k = 0,l,.... 

In view of Lemma 13.31 a natural idea is to relax the "difficult" constraint 
1 -9 e C 2d (K) in ([31]) to M h (l-g, z) ^ for fixed k, and then let fc oo. 
Indeed, for every fixed k, the latter is much easier to handle as it defines 
a spectrahedroro C M^ 2d ) on the coefficients of the homogeneous poly- 
nomial g. And so one obtains a hierarchy of convex relaxations of V by 
minimizing the (strictly) convex function g i— )■ J exp(-g) on the spectrahe- 
dra Afc, /c G N, which yields a monotone nondecreasing sequence of lower 
bounds pk < p, k G N, on p. Of course the larger k (i.e., the larger the size 
of the localizing matrix Mfc(l — g, z)) the better the lower bound pk (but 
also the harder the problem). 

Theorem 3.4. Let K C f" be compact with nonempty interior and consider 
the finite- dimensional convex optimization problem: 

(3.4) V k : p k = inf ( / exp(- 5 )dx : M fc (l- 5 , z) b 

9eP[x] 2 d I JR™ 

T/te sequence (pt), k G N, is monotone nondecreasing with pk < p and: 

(a) // Pfc /ias an optimal solution g* k G P[x] 2 d i/ien £/iere exists a SOS 
polynomial a* k G £[x]fc such that 

(3.5) f x a exp(-g* k )dx = [ x a a* k dp, V|a| = 2d. 

(3.6) / (l-g* k )atdv = 0, 

and p k =™ cr* k dp. 
JVL 

(b) Conversely if h3. 5\) - [37S\) holds for some g k G P[x] 2 ^ and some a k G 
S[x]fc then g k is an optimal solution ofVk- 

(c) If in addition, 

sup sup \g* ka \ < M, 

k \a\=2d 

for some M > 0, then pk —> p and V has an optimal solution g* . 

For a proof see §4.11 So problem Vk amounts to minimize a strictly convex 
function on a spectrahedron of M^ 2a! ) . For instance, one may use interior 
point methods and minimize the standard log-barrier function 

g^4>u{g)-= ( exp(-g)dx- ~logdet(M fe (l - g, z)), 

JR™ v 

spectrahedron is a convex set that can be formed by intersecting the cone of positive 
semidefinite matrices S n with a linear affine subspace. 
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with parameter v, and let v — > oo. For more details on log-barrier methods, 
the reader is referred to e.g. Wright [25J. 

3.2. Upper bounds via outer approximations. Let K C W 1 be the 

compact basic semi-algebraic set defined by 

K = {x : itj(x) > 0, j = l,...,m}, 

for some polynomials (v,j) C M[x]. As K is compact, with no loss of gen- 
erality we may and will suppose that ui(x) = M — ||x|| 2 for some M large 
enough. Let 

m 

Q* = {^2<JjUj : uj E S[x], j = 0, . . . ,m}, 

3=0 

be the quadratic module of R[x] generated by the Uj's (with uq = 1). Q* is 
Archimedean because the quadratic polynomial x i— > M — ||x|| 2 belongs to 
Q* . If one defines 

m 

Q* k := {^2&jUj : Uj £ S[x]; degajUj < 2k, j = 0, ... , m}, k € N, 

3=0 

(a subset of Q* with a degree bound on the SOS weights <7j) then this time 
one may replace the "difficult" constraint 1 — g £ C2d(K) by the stronger 
constraint 1 — g G Q^. for fixed E N. 

The convex cone Q* k is the dual cone of the closed pointed convex cone 

(3.7) Q k := {y E M s(2fc) : M fc (y) ^ 0, M fc _„.(^-, y) ^ 0, j = 1, . . . ,m}, 

where Mfc(y) (resp. Mfc_„. (uj, y)) is the moment matrix associated with 
the sequence y (resp. the localizing matrix associated with y and the poly- 
nomial Uj). Hence, Q* k has a nonempty interior; see e.g. Rockafellar [21J. 
Then solving the problem 

(3.8) V' k : inf { / exp(- 5 )dx : 1 - g E Q% ) , k E N, 

k> d, now provides a monotone sequence of upper bounds p' k > p, k E N. 

Theorem 3.5. Let K C R n 6e compact with nonempty interior and consider 
the finite- dimensional convex optimization problem t3.8\) . k > d. 

(a) The sequence (p' k ), d < k E N ; is monotone nonincreasing with p' k — > p 
as k —J- oo. Moreover there exists an optimal solution gt E P[x]2d. 

In addition, assume that there exists go E P[x]2d such that l — go£ int Q k 
for some k$ > d. Then: 

(b) If k > ko and g k E P[x]2^ is an optimal solution ofVi, there exists a 
vector y k E Qk such that: 

(3.9) 0=(l-^,y fc ); y k a = f x a exp(-^) dx, V \a\ = 2d. 
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(c) Conversely, if k > and (y fc , 1 — g k ) £ Qk x satisfy A3.S\) . then 
g* k is an optimal solution of V' k . 

(d) Let (gl) C P[x]2d, k > ko, be a sequence of optimal solutions. If 

sup sup \g* ka \ < N, 

k \a\=2d 

for some N > 0, then every accumulation point g* of the sequence (gt)> 
k € N, is an optimal solution ofV. 

For proof see §4.21 

4. Proofs 

4.1. Proof of Theorem [3741 

Proof, (a) That the sequence (pk) is monotone non decreasing is straightfor- 
ward since the constraints of Vk are more and more restrictive as k increases. 
Next, as K is compact there exists M,5 > such that M — ||x|| M > 5 for 
all x € K. With go £ £[ x ]<2 being the polynomial x i— > Af _1 ||x|| 2c( , one has 
1 — 5o( x ) > S/M for all x G K and so M&(1 — go, z) y 0. Indeed, for every 
f € R[x] fc , 

(f,M fc (l- go ,z)f) = J f(l-g )dp > A J^/ 2 ^ > 0. 

where the last inequality is because K has nonempty interior and supp p = 
K. Observe also that go G Pfxkd and so go is a strictly feasible solution of 
Vki that is, Slater's conditiory holds for Vk- So the Karush-Kuhn- Tucker 
(KKT) optimality conditions at a point g k are necessary and sufficient for 
g k to be a (global) minimizer of Vk- Therefore, there exists 0^ As <S s (k) 
such that: 

(4.1) f -x a exp(-^)dx + (A,M fc (x Q , z)> = 0, V|a| = 2d, 
and 

(4.2) (M d (l -g* k , z),A) = 0, 

where one has written Mfc(<?^ z) = 9% a Mfc(x a , z) (with Mfc(x a , z) being 
the localizing matrix with respect to the polynomial x a and the moment 
sequence z). From its spectral decomposition, A = ^ ■ qjqj for some vectors 
(c/j), which yields 

(A,M fc (x Q , z)) = I aU*)x a dp, V|a| = 2d, 
Jk 



For a convex optimization inf x {/(x) : ht(x) > 0, k — 1, ...,m} Slater's condition 
holds if there exists xo such that ftfe(xo) > for all k. In this case the Karush-Kuhn- Tucker 
(KKT) optimality conditions are necessary and sufficient. 
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where a* k = ]Cj(</J v fc( x )) 2 G is a SOS polynomial. And so (|4.1h yields 
(|3.5p . and the complementary condition (|4.2|) yields (|3,6|) . Next, multiplying 
both sides of (|4.ip with ^ a and summing up yields: 

/ 9k exp(-fffc)dx ( = — / exp(-^)dx) = (A, M k (g* k , z)) 

a* k g* k dp, 

K 



/ 4 dp [by||OJ) 



(b) The converse is because under Slater's condition, the KKT optimality 
conditions are also sufficient for gt to be a minimizer. 

(c) Write g* k (x) = T,\ a \=2d 9* ka x a ■ As sup fc sup H=2(i < M, there 
exists a subsequence (fcj), i e N, and a vector g* G such that for every 
a G N n with |a| = 2d, g* k _ a -> 5* as i -> 00. 

Notice that x h-> <?*(x) = limsup^oo ^.(x) = liminf g^. (x). Next, as 
P > Pk for every fc, and exp(— g k ) > 0, by Fatou's Lemma (see e.g. Ash |2j) 



p > liminf p k . = liminf / exp(— gt.) (fx 
> / liminf exp(—<7^ ) dx 

= / exp(-limsuppfc.)dx = / exp(-£i*)dx. 

On the other hand, observe that for every k G N, M/u(l — gj, z) ^ implies 
rvlj(l — 5^, z) ^ for all j < A;. Hence, let j be fixed arbitrary so that 
Mj(l — z) ^ 0, for all sufficiently large k. For every / G M[x]j, f 2 g k is 
uniformly bounded on K and so by Fatou's Lemma 

0<limsup/ f 2 (l-gljdp = -liminf / f 2 (g* ki - 1) dp 

< / -liminf/ 2 (g* k . - I) dp 
Jk 



[ f 2 (l-g*)dp. 
JK 



As / G M[x] j was arbitrary, Mj(l — g* , z) ^ 0; and as j was arbitrary, 
Mj(l — g* , z) ^ 0, for all j = 0,1,... But by Lemma [3.31 this implies 
1 — g* > on K, and so g* is a feasible solution for V . Combining with 
p > J" Rn exp(— g*)dx yields that g* is an optimal solution of V . □ 

4.2. Proof of Theorem \3M 

Proof, (a) The sequence (p' k ) is monotone non increasing because Q* k C Q* k , x 
for every A; G N. Next, let e > be fixed, and let G P[x] 2 d be such that 
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P < /gn exp(— g)dx < p + e. Observe that g := (1 + e) g G P[x]2d and 
/ exp(-g)dx = / exp(-(l + e)~ l g) dx 

= (1 + e) n ' d [ exp(- 5 ) dx < (1 + e) n / d (p + e). 

JR™ 

and 1 — g = > on K. Therefore, as Q* is Archimedean, by Putinar's 

Positivstellensatz [2D], 1 — g G Q*, that is, 1 — g G Q| for some k\. Hence 
g is a feasible solution of Vk x which implies p < p' k < (1 + e) n / d (p + e). As 
the sequence is monotone and e > was arbitrary, we obtain the desired 
convergence p' k — > p. 

Next, let w = (w a ), a G N n , be the sequence of moment of the Lebesgue 
measure on K, and let g G R[x]2d be an arbitrary feasible solution. Let 
Mfc(w) (resp. Mfc_ v .(uj w)) be the moment (resp. localizing) matrix associ- 
ated with w and Uj. Hence, M^(w) >- and M^_ v . {v,j w) >- 0, j = 1, . . . , m, 
because K has nonempty interior. The constraint 1 — gG Qk( u ) which reads 
(1 — g) = cr + ]Cjli ^""j f° r some SOS polynomials (<Tj) has the equivalent 
form 

m 

(4.3) (l-^)a = (X ,B a )+^(X i ,C ia ), Va€NJfc, 

i=i 

for some appropriate real symmetric s(2(k — Vj)) x s(2(k — Vj)) matrices 
Xj y 0, (and where B a , Cj Q are given real symmetric matrices). The vec- 
tor of coefficients of the SOS polynomial aj G S[x] is obtained from the 
eigenvectors of Xj, j = 0, . . . , m. For more details see e.g. [11] . Multiplying 
(14. 3p by w a and summing up yield: 

„ m 

(l-g)dx = (X ,^2~B a w a ) + ^2(Xj,^2c ja w a ) 

^ K a j=l a 

m 

= (X ,Mfc(w))+5^(X ij Mfc_ 1)J (u i w)) 

3=1 

Observe that every feasible solution g G M[x]2d is nonnegative otherwise 
j exp(—g)dx is not bounded. And so exp(— g) > 1 — g on K because g > 0. 
Therefore, for any minimizing sequence (g n ,X") of (|3.8p . 



/ (l-g n )dx< / exp(-# n )dx < / exp(-g n )dx < / exp(-sr )dx. 

JK JK JR" JR" 



And so 
sup < 



(Xg,M fc (w))+y;(Xy,M fc _ l>j (« i w)) I < / exp(-g )dx. 



As Mfc(w), Mk-Vj ( u j w ) >~ and X" ^ 0, j = 1, . . . , m, all matrices X" are 
uniformly bounded. Hence one may extract a subsequence {rig) such that 
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X™ £ — >• Xj; >z 0, j = 0, . . . , n, as I — > oo. And so for every a G N^., 

(X*,B a ) + V(X),C iQ ) = lim (1 =: (1 -<?*%, 

for some homogeneous polynomial g* G M[x]2d (as all coefficients g* a with 
|q| ^ 2<i vanish). In other words {fjng)a 

->■ g* as £ — >■ oo, for all a G N^,. 
As (1 — g ni ) > on K one also has (1 — g*) > on K. Finally, since we 
also have the pointwise convergence g ne ( x ) 5*( x ) f° r an x £ invoking 
Fatou's lemma yields 



Pk 



lim / exp(— g n( ,) <ix > / liminf exp(— g ne ) dx 

'^°°JR™ ' ' ! JR™ 

= / exp(-5*)dx, 

JR™ 



which proves that g* is an optimal solution of (|3,8p . 

(b) 1 — <7o G m t (5^ for all k > kg because 1 — go G hit Qjjl o and Q* k D Ql . 
Hence Slater's condition holds for "Pi whenever k > ko. Therefore, if g k is an 
optimal solution of Vk, there exists y k G Qk such that the KKT-optimality 
condition hold, which yields (|3.9p . 

(c) Follows from the fact that under Slater's condition, the KKT-optimality 
conditions are sufficient for g* k to be an optimal solution of V' k . 

(d) If sup fc sup| Q ,| =2 d \g ka | < N, let g* be an accumulation point, i.e., 
a limit point of some subsequence (g k ), i G N, i.e., such that for every a 
with | a | = 2d, g k . a . — > <7* as i — > oo. Let g* G M[x]2d be the homogeneous 
polynomial with coefficients g* a , |a| = 2d. For every x G K, (1 — ^(x)) > 
because 1 — <^ . G Q* k . for every i. Therefore, with x G K fixed, 

0< 1- lim <£(x) = l-ff*(x), 

I— s>oo 

and so, as x G K is arbitrary, 1 — g* G C2<f(K), which implies that 1 — 5* is 
a feasible solution of P. Next, from (a) and by Fatou's Lemma, 

p = lim p ki = liminf^ = liminf / ex.p(-g* k )dx 

> / liminf exp(— g£.) cfoc 
= / exp(-£*)dx, 

JR" 

and so as 5* is feasible for "P, it is an optimal solution of "P. □ 

4.3. The dual cone of Cd(K). Recall that M(K) is the space of finite 
Borel measures on K. Let Coo(K) C M[x] be the space of polynomials 
nonnegative on K C M. n and let Aoo C M[x]* be the set: 

(4.4) := j (J x Q dcf^J , a G N n : <fi G M(K 
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It is known that when K C W 1 is closed then A^ = Coo(K) and Aoo = 
Coo(K)*. See e.g. the proof in Laurent [161 Prop. 4.5] for K = R n , which 
also works for any closed set K C M. n . The following result is a truncated 
version when K is compact. 

Lemma 4.1. Let K C M. n be compact. For every d G N, the dual Crf(K)* 
o/C,i(K) is the set: 

(4.5) A d := |QV#),aGN2 : <^>G M(K) j . 

Proof. For every y = G A^ and / G C d (K) with coefficient vector 

f G 

(4.6) (y,f) = f*y a = I = I f d( t>> 0- 

Since f|4.6|) holds for all / G C d (K) and y G A^, then necessarily A^ C 
C d (K)* and similarly, C d (K) C AJ. Next, 

A^ = {f el s M : (f,y) >0 VyG A,} 

/€R[x] d : j^/#>0 V0GM(K) 
A^CQ(K), 

and so A d = C d (K). Hence the result follows if one proves that A^ is closed, 
because then C d (K)* = (A* d )* = A d , the desired result. So let (y fc ) C A d , 
k G N, with y k — > y as k — > oo. Equivalently, f K x a d(j)k — > y a for all a G N^. 
In particular, the convergence y$ — > yo implies that the sequence of measures 
{4>k)i k G N, is bounded, that is, sup fc <^fc(K) < M for some M > 0. As K 
is compact, the unit ball of M(K) is sequentially compact in the weak * 
topology <r(M(K), C(K)) where C(K) is the space of continuous functions 
on K. Hence there is a finite Borel measure 4> £ M(K) and a subsequence 
(ki) such that j* K gci^ — > / K ffd^ as i — > oo, for all g G C(K). In particular, 
for every a G NJJ, 



y Q = lim = lim y% = lim / x° 



X 
K 



which shows that y G A^, and so A^ is closed. □ 
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